Short-time dynamic in the Majority vote modeh The ordered and disordered initial 

cases 
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This work presents short-time Monte Carlo simulations for the two dimensional Majority-vote 
model starting from ordered and disordered states. It has been found that there are two pseudo- 
critical points, each one within the error-bar range of previous reported values performed using 
fourth order cumulant crossing method. The results show that the short-time dynamic for this 
model has a dependence on the initial conditions. Based on this dependence a method is proposed 
for the evaluation of the pseudo critical points and the extraction of the dynamical critical exponent 
z and the static critical exponent P/v for this model. 
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I. INTRODUCTION 



Critical phenomena in equilibrium statistical systems 
is one of the most important topics in physics. Much 
of the attention has been focused on the universality of 
the critical exponents, with several universality classes al- 
ready characterized in equilibrium systems. On the other 
hand the critical behavior of non-equilibrium statistical 
systems has been receiving a lot of attention in recent 
years, but the characterization of the different universal- 
ity classes is far from be complete. One of the simplest 
non-equilibrium models is the two dimensional majority 
vote model, an Ising-like system (up-down symmetry and 
spin-flip dynamic) with a continuous order-disorder tran- 
sition with the same critical exponents that the two di- 
mensional Ising modelJlHsl , as expected from the predic- 
tion of Grinstein aZ [4| : every spin system with spin-flip 
dynamic and up-down symmetry falls in the Ising model 
universality class. However, there is some controversy 
about the universality class for higher dimensions. A re- 
cent work claims Q that the upper critical dimension for 
the majority vote model is 6 instead of 4, based on nu- 
merical calculations. Another discrepancy in the critical 
exponents have been found in simulations on non-regular 
lattices [1,0]- It must be mentioned that all of the results 
mentioned above were performed using standard " Monte 
Carlo" simulations and Finite Size Scaling approaches 
for the evaluation of the static critical exponents. On 
the other hand the time evolution can gives important 
information about the universality of a given system. It 
has been shown by Janssen et al. [1] that when systems 
with relaxation dynamics are quenched from high tem- 
peratures to the critical temperature there is a short crit- 
ical universal behavior. Numerical simulations have con- 
firmed this behavior in the Ising and the Pott models (see 



reference [ol for a review of these results). Concerning 
the critical dynamic in systems without detailed balance, 
there are some works that evaluate the critical dynamic 
exponent z Bead the fluctuation-dissipation ra- 

tio Xqo for the Majority model, but always starting 
from a disordered state. As expected the results were 
compatibles with the Ising ones. Given all these results 
one should expect that the basic assumption for the dy- 
namic relaxation of the k-th moment of the order param- 
eter starting from an completely ordered state will be the 
same as in the Ising model, but this has not been proved 
yet. 

The aims of this work are: a) to evaluate the critical 
point using short time dynamic starting from ordered and 
disordered initial conditions and b) evaluate the dynamic 
critical exponent z and the static critical exponent (3 /v. 
This will test if the Grinstein prediction holds for the 
short time dynamic in the majority vote model. 



II. MODELS AND DEFINITIONS 

In the Majority vote model each lattice site has a spin 
whose values are tr = ±1 and its dynamic can be grouped 
by the spin flip rule 



W, = \[l-a,,f{H,)], 



(1) 



here Hi is the local flcld X^nn '''i 07=^2, ±4 produced 
by its nearest neighbors and / is a function with up 
down symmetry that depends on two control parameters 
/(O) = 0, /(2) = -/(-2) = X and /(4) = -/(4) = y. 
The parameters x and y can be associated with inter- 
face and bulk temperatures respectively [iQ] using the 
relations 



tanh2/32, j/ = tanh4/34. 



(2) 
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The Majority model is obtained setting x — y {(3^ < 
/32)and the critical point is at Xc = 0.850(2) §]. 



2 



The equilibrium case can be obtained along the line 
y — 2x/(l + x^) (Glauber dynamic), where the temper- 
atures are equal (/32 — Pa — P) and the critical point is 
= ilog(l + V2). 

The order parameter is the standard for Ising-like sys- 
tems, defined by 

"i=^(Xl'^») (3) 

i 

where N = \s the total number of lattice sites [L is 
the lateral size. 

Starting from a disordered state the dynamic for the k- 
th moment of the order parameter was derived by Janssen 
the mathematical expression is given by 

(4) 

here t = [T — Tc)/Tc is the reduced temperature, or the 
reduced control parameter for non-equilibrium systems, 
t is the dynamic time variable, b is the re-scaling factor, 
(5 and v are static critical exponents, z is the dynamic 
critical exponent, is the scaling dimension of the ini- 
tial (small) order parameter toq (see and reference 
therein) . 

At the critical point for sufficiently small mp and large 
systems (L — >■ oo) the order parameter follows a power 
law dynamic 

m{t) - mo^^ (5) 

where 9 is defined by {xq — j3/v)/z. One must remark 
that there is a strong dependence on the initial value of 
the order parameter. This dependence does not affect the 
power law behavior at the critical point, what is affected 
is the exponent 9 that tends to the real value in the limit 
mo — >■ 

For the dynamic starting, from the ordered state we 
have the assumption that the scaling dynamic form is 
given by 

m^''\t,T,L) = b'^'^/''m'^''\b-Hy/''T,b-~^L), (6) 

again it can be obtained the scaling form at the critical 
point taking the limit L oo. 

M{t) - t-^/"''. (7) 

This has been proved in equilibrium systems, like the 
Ising or the Potts model. 

For the evaluation of the critical point it has been used 
the fact that theoretically the order parameter evolves 
as a power law at the critical point If it is evaluated the 
difference between a power law and the time evolution for 
different values of the control parameter x we will expect 
a minimum in those differences at the critical point. 

The simulations were performed choosing the lattice 
site randomly starting from both, a completely ordered 
state (m = 1.0) and a carefully prepared disordered state 
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FIG. 1: (color online) a) Dynamic starting from an ordered 
state, b) Dynamic starting from a disordered state (mo = 
0.0875). In both cases the top continuous line corresponds to 
a X above the critical point and the lower one to a a; below 
the critical point. The dashed line shows the expected power 
law behavior. 



with small magnetization values toq < 0.1. The or- 
der parameter is evaluated as a function of the time 
with At = 1 corresponding to a Monte Carlo time step 
(MCTS). In order to avoid finite size effects we used lat- 
tices with lateral size L — 2^ [N — L?). The average were 
taken with at least 10^ independent simulations and 250 
and 1000 MCTS for the evaluation of the critical point 
and dynamical exponents respectively. 



III. RESULTS 

Given that the power law behavior of the order pa- 
rameter as function of t (decay from an ordered state, 
increase for a disordered one) one can performed simula- 
tions for different values of x around the expected critical 
point. In Fig. [T] we can observe that above and below a 
certain value of x we have dynamic that differ from a 
power law (illustrated by the dashed line). 

From here one must define the criterion that will mea- 
sure which is the best power law curve, in this work it 
was used the measurement of the for each curve with 
respect to a power law behavior. I must remark that pre- 
vious results for equilibrium systems give the same result 
for the critical point for both initial states. However, two 
different values were found for the majority vote model, 
see Fig. H] 

The critical point values obtained were Xc = 
0.85007(6) and Xc = 0.85147(2) for the order and disorder 
state respectively, both results are in perfect agreement 
with the obtained in the static case (references |lH3|) and 
are quite close. 

We have a clearly dependence on the initial condition 
in short time dynamic of this model, which is not the 
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FIG. 2: (color online) Deviation from the power line behavior 
for the decay (mo = 1, left curve) and the increase (mo = 
0.0875, right curve) CEises. 



TABLE I: Critical point for initial disordered states. 
mo I 0.0375 0.0500 0.0625 0.0750 0.0875 
Xc 0.84929(10) 0.84972(9) 0.85019(6) 0.85076(5) 0.85147(2) 



case for equilibrium systems. There is no doubt about 
the critical point obtained from the ordered phase, since 
Too = 1 is one of the fixed point under renormalization 
group transformations. However, one must remember 
that the power law is valid only in the limit toq — >■ (the 
other fixed point) for the disorder case, assuming that 
the "real" critical point is located at this limit, one can 
proceed to evaluate the critical point for another values 
of Too and with this values extrapolate the critical point 
for the disordered phase. The results are showed in tabltU 



Once that each value Xc(too) has been evaluated, the 
dynamical exponent Q can be obtained. For the evalua- 
tion of this exponent the simulations were performed with 
1000 MCTS discarding the first time steps, since there is 
an initial time scale tmic where the power law stabilizes 
triiic ~ 20 (see Fig. [3]). The results for the exponent d are 
showed in table HIl 



With a extrapolation of these values to too = the 
value of the critical point and the d exponent were eval- 
uated (see Fig. U]). The result for the critical point was 
Xc = 0.84860(10), which is clearly different from the or- 
dered one, however both values are within the error bar 
from the obtained in the static case 0.848 < Xc< 0.852. 
From now on the pseudo critical point evaluated with 



TABLE IL theta exponents for growing process, 
mo 0.0375 0.0500 0.0625 0.0750 0.0875 




FIG. 3: (color online) Growing of the order param- 
eter at a;c(mo), the continuous curves are for mo = 
0.0875, 0.075, 0.0625, 0.5, 0.0375 from top to bottom. The 
dashed line represents the power law growing with 9 — 0.1751 
(the result in this work for the majority vote model). 



the decay process will be denoted as x° and the evalu- 
ated with the growing process with x^. This surprising 
result seems similar to the obtained for weak first-order 
phase transitions [lO^, where two pseudo-critical points 
exits due to the metastable states above and below the 
critical point. However, there is an important difference 
in this case: for weak order phase transitions the smaller 
critical point corresponds to the decay process, and the 
bigger corresponds to the growing process, contrary to 
the majority vote model case. 

In the evaluation of the exponent 6 a linear extrapola- 
tion gives the result of 0.175(3), which is lower from the 
values of the two dimensional Ising model, 6 = 0.191(1) 

ni, 



and from the previously evaluated in references 
for the majority vote model, 9 = 0.192(2). The difference 
with respect to the Ising model could be understood con- 
sidering that the results obtained here seems to indicate 
a hole new dynamic. The differences with previous re- 
sults for the majority model can be explained observing 
the simulations details used previously: first the critical 
point used was a;=0.850, which is above the result for 
x'^. Second the systems sizes used previously were really 
small (L — 32), at this size the growing process is not 
very long and is really hard to see the power law behav- 
ior. 

One can obtain the dynamical exponent z evaluating 
the second moment of the magnetization at the critical 
point x^ 



,(2) 



ty, y = {d-2l3/v)/z, 



and the autocorrelation 



e 0.1769(8) 0.1774(4) 0.1782(3) 0.1788(4) 0.1792(3) 



(8) 



(9) 
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a) 



b) 



TABLE III: Summary of the results in this work and of the 
Ising model. 



FIG. 4: a) Evaluation of the critical point Xc and b) the 
exponent. 
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FIG. 5: (color online) a) Evaluation of A, the continuous line 
shows the autocorrelation time evolution and the dashed line 
shows the power law behavior with A = 0.758. b) Evaluation 
of y , the continuous line shows the second moment order 
parameter and the dashed line shows the power law behavior 
with y = 0.799. 



Both startmg from mo — and using 1000 MCTS. Again 
there is a tmic in each case (around 20 for the autocor- 
relation and 75 for the second moment, see Fig. [S]). The 
results obtained are y = 0.799(17) and A = 0.758(2). 
Combining both results it can be obtained the values 
z = 2.143(9) and P/f = 0.143(18). A summary of the 
results are showed in table IIIIl where it can be observed 
discrepancies between most of the values for the major- 
ity model and the Ising ones. It must be remark that all 
these results were obtained using just the growing pro- 
cess. 



Finally the decay exponent fi/i'z was evaluated start- 
ing from an ordered phase (at x°), using 1000 MCTS 
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FIG. 6: (color online) Order parameter relaxation (mo = 1), 
here the dashed line shows the power law behavior with P/v = 
0.0526. 



(Fig. inj, the result was ji/v = 0.0526(5), that is lower 
compared to the Ising one, 0.0580(5). Again the first 
time steps were discarded for the evaluation of the expo- 
nent (Fig. [B]) . Theoretically it is possible to obtain the z 
exponent using the known value of or knowing the z 
value one can obtain the P /v value, but in both cases the 
results depend on values obtained with a growing process 
(z) or the static simulations {/3/nu). The approach taken 
in this work is that the growing and the decay process 
are different and it could be possible that the dynamic 
exponent z is different in each case, so for the decay case 
the reporting value is z = 2.37(2). 

The fact that we have two pseudo-critical points (that 
not corresponds to a weak phase transition) and that the 
decay and growing process are slower that in the Ising 
model must be related to the absence of detailed balance 
condition. One of the consequences of this absence is 
that we do not have a unique thermodynamic temper- 
ature, in this case we have two, so looking at the snap- 
shots for different initial conditions at the pseudo-critical 
points we can speculate about the competition between 
the two "temperatures" that governs the dynamic in non- 
equilibrium Ising systems. Figure [7] shows the time evo- 
lution with Too = at the two pseudo-critical points, a) 
for x'^ and b) for x°. Initially the number of sites with 
spin-flip probability depending on /32 are very similar to 
the ones depending on /?4, the time increases from left to 
right and it can be observed that at x° the coarsening 
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FIG. 8: Snapshots at the pseudo-critical points starting with 
mo = 1, a) and b) x°. Times are (from left to right) 1, 
100, 1000, 10000 and 20000. 

seems to appear faster that at . 

Figure |S] shows the same time evolution for mo = 1. 
In this case at the beginning of the evolution all sites 
have spin-flip probability that depends only on (3/^ and 
the decay is slightly lower at x°. 

It seems that the coarsening differs at the two pseudo- 
critical points, so as a final test the time evolution for 
a special initial condition were performed setting mo re- 
ally close to zero putting almost half of the spins in one 
state and the other half in the other state with a circular 
border, in this way we have a large number of sites with 
/34 while the number of possible sites with 62 increases 
(Fig. in]). We can observe that the circular shape last 
longer at x°. In order to corroborate the effect of the 
difference between temperatures it should be performed 
simulations in spin-like systems for different ratios (34/(32, 
except for the equilibrium case (^4/(^2 = 1- 



IV. CONCLUSION 



In this work it has been shown that the short time 
dynamic in the majority vote model presents power law 
behavior at different control parameters for the grow- 
ing, — 0.84860(10), and the decay processes, x° = 
0.85007(6). These pseudo-critical points are compati- 
bles with results for the critical point reported previously, 
Xc — 0.850(2). It has been show also that the dynamic 
in both cases is slower that in Ising model for all the 
quantities calculated (m, m^^^ and A). These results 
seems to be related to the competing dynamic between 
the interface (/32) and bulk {(34) temperatures associated 
to the dynamic, and as consequence to the absence of 
detailed balance in the system. In order to corroborate 
these results additional simulations must be carry on in 
systems without detailed balance. The dynamical critical 
exponent (z) and the static critical exponent {(3/i^) has 
been evaluated independently using a growing process, in 
both cases the results were close to the Ising ones. For 
the decay process the z exponent was evaluated using re- 
sults from static simulations founding that the value is 
different from the obtained in the growing process. 
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FIG. 9: Snapshots at the pseudo-critical points starting with 
mo almost zero and a circular border, a) and b) x°. Times 
are (from left to right) 1, 100, 1000, 10000 and 20000. 
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